Quantum phase transitions in the sub-ohmic spin-boson model: 
Failure of the quantum— classical mapping 
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The effective theories for many quantum phase transitions can be mapped onto those of classical 
transitions. Here we show that the naive mapping fails for the sub-ohmic spin-boson model which 
describes a two-level system coupled to a bosonic bath with power-law spectral density, J {to) oc ui s . 
Using an e expansion we prove that this model has a quantum transition controlled by an interacting 
fixed point at small s, and support this by numerical calculations. In contrast, the corresponding 
classical long-range Ising model is known to display mean-field transition behavior for < s < 1/2, 
controlled by a non-interacting fixed point. The failure of the quantum-classical mapping is argued 
to arise from the long-ranged interaction in imaginary time in the quantum model. 



Low-energy theories for certain classes of quantum 
phase transitions in clean systems with d spatial dimen- 
sions are known to be equivalent to the ones of classi- 
cal phase transitions in (d + z) dimensions, where z ist 
the dynamical exponent of the quantum transition [1]. 
This mapping is usually established in a path integral 
formulation of the effective action for the order parame- 
ter, where imaginary time in the quantum problem takes 
the role of z additional space dimensions in the classical 
counterpart. The tuning parameter for the phase transi- 
tion, being the ratio of certain coupling constants in the 
quantum problem (where T is fixed to zero), becomes 
temperature for the classical transition. For the quan- 
tum Ising model, where the transverse field can drive the 
system into a disordered phase at T = 0, the quantum- 
classical equivalence in the scaling limit can be explicitly 
shown using transfer matrix techniques [1]. While this 
formal proof is only applicable for short-range interac- 
tions in time direction, it is believed that it also holds for 
long-range interactions, which can arise upon integrat- 
ing out gapless degrees of freedom coupled to the order 
parameter. (Counter-examples are phase transitions in 
itinerant magnets, where the elimination of low-energy 
fermions produces non-analyticities in the resulting or- 
der parameter field theory [2].) A paradigmatic example 
is the spin-boson model [3, 4], where an Ising spin (i.e. 
a generic two-level system) is coupled to a bath of har- 
monic oscillators: eliminating the bath variables leads 
to a retarded self-interaction for the local spin degree of 
freedom, which decays as 1/t 2 in the well-studied case 
of ohmic damping. Interestingly, the same model is ob- 
tained as the low-energy limit of the anisotropic Kondo 
model which describes a spin- 1/2 magnetic impurity cou- 
pled to a gas of conduction electrons [5, 6]. 

The purpose of this paper is to point out that the naive 
quantum-classical mapping can fail for long-ranged inter- 
actions in imaginary time even for the simplest case of 
(0 + 1) dimensions and Ising symmetry. We shall explic- 
itly prove this failure for the sub-ohmic spin-boson model, 



by showing that the phase transitions in the quantum 
problem and in the corresponding classical long-range 
Ising model fall in different universality classes. 

The spin-boson model is described by the Hamiltonian 



H S b = -^Px+^v z +y^ j ijJia\ai+ C ^-y^ j Xi(ai+a\) (1) 



in standard notation. The coupling between spin a and 
the bosonic bath with oscillators {a^} is completely spec- 
ified by the bath spectral function 



J (w) = 7T ^ ( W _ ^i) > 



conveniently parametrized as 

J(w) = 2na ujI~ s oj s , < to < oj c 



(2) 



s > -1 (3) 



where the dimensionless parameter a characterizes the 
dissipation strength, and oj c is a cutoff energy. The value 
s = 1 represents the case of ohmic dissipation, where 
a Kosterlitz-Thouless transition separates a delocalized 
phase at small a from a localized phase at large a. These 
two phases asymptotically correspond to eigenstates of 
o~ x and cr z , respectively. 

In the following, we are interested in sub-ohmic damp- 
ing, < s < 1 [12, 13]. The standard approach is to 
integrate out the bath, leading to an effective interaction 



<5int. — 



J dTdT'<J z (T) 9 (T-T')a z (T') 



(4) 



with g(r) oc 1/t 1+s at long times. Numerical renormal- 
ization group (NRG) calculations in Refs. 7, 8, performed 
directly for the sub-ohmic spin-boson model, have estab- 
lished that a second-order quantum transition occurs for 
all < s < 1. Here we use an analytical renormalization 
group (RG) expansion, controlled by the small parameter 
s, to establish that the spin-boson transition at small s is 
governed by an interacting fixed point with strong hyper- 
scaling properties. This analytical result is supported by 
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NRG calculations. In contrast, the transition in the clas- 
sical Ising model is known to display mean-field behavior 
for < s < 1/2 [9, 10]. 

Scaling and critical exponents. A scaling ansatz for the 
impurity part of the free energy takes the form 

F imp = Tf{\a-a c \T- 1 '»,eT-y') (5) 

where \a — a c \ measures the distance to criticality. The 
bias e takes the role of a local field (with scaling exponent 
y e ); and v is the correlation length exponent which de- 
scribes the vanishing of the energy scale T* , above which 
critical behavior is observed [1]: T* cx \a — a c \ v . The 
ansatz (5) assumes the fixed point to be interacting; for 
a Gaussian fixed point the scaling function will also de- 
pend upon dangerously irrelevant variables. 

With the local magnetization M\ oc = (a z ) = 
—dFi mp /de and the susceptibility xioc = —d 2 Fi mp /(de) 2 
we can define critical exponents (see also Ref. 16): 

Mi oc (a> a c ,T = 0,e = 0) cx (a-a c f, 

Xioc(a < ct c ,T = 0) cx (a c - a)~ 7 , 

M loc (a = a c ,T = 0) cx | e | 1 / ^ , 

Xioc(a = a c ,T) cx T~ x , 

x" oc (a = a c ,T = 0,uj) cx \u\~ v sga(u). 

The last equation describes the dynamical scaling of xioc- 
In the absence of a dangerously irrelevant variable 
there are only two independent exponents, e.g., v and 
y e . The scaling form (5) yields hyperscaling relations: 

/3 = 7^, 2/3 + 7 = ^, 7 = ra , J=i±|. (7) 

Hyperscaling also implies x — y, which is equivalent to 
so-called lo/T scaling in the dynamical behavior. 

Long-range Ising model. The classical counterpart of 
the spin-boson model (1) is the one-dimensional Ising 
model [3, 4] 

Hcl = - Ji]SiSj + ^sr (8) 

m 

with interaction = J/\i — j\ 1+s . Hsr contains an ad- 
ditional generic short-range interaction which arises from 
the transverse field, but is believed to be irrelevant for 
the critical behavior [9, 10]. As proven by Dyson [14] this 
model displays a phase transition for < s < 1. Both an- 
alytical arguments, based on the equivalence to a 0(1) </> 4 
theory [9], and extensive numerical simulations [10] show 
that the upper-critical dimension for the d-dimensional 
long-range Ising model is c?+ = 2s, i.e., in d = 1 the tran- 
sition obeys non-trivial critical behavior for 1/2 < s < 1. 
In contrast, mean-field behavior obtains for < s < 1/2, 
with exponents [9, 10] (3 = 1/2, 7 = 1, 6 = 3, v = 1/s, 
y = s violating hyperscaling. 



As the exponent s exclusively determines the power 
laws of spectra and correlations in a long-range model 
once spatial dimensionality (d = 1) is fixed, s takes the 
role of a "dimension", i.e., we will refer to s = 1/2 as 
upper-critical "dimension" of the classical Ising chain. 

Near s = 1 the phase transition can be analyzed us- 
ing a kink gas representation of the partition function, 
where the kinks represent Ising domain walls [15]. This 
expansion, controlled by the smallncss of the kink fugac- 
ity, is done around the ordered phase of the Ising model, 
corresponding to the localized fixed point of the spin- 
boson quantum problem. For small (1 — s) the results 
obtained via the perturbative kink-gas RG are consis- 
tent with the NRG data for the spin-boson transition [7] , 
indicating that the quantum-classical mapping works in 
the asymptotic vicinity of the localized fixed point. 

Spin-boson model: perturbative RG. We now describe 
a novel RG expansion which is performed around the 
delocalized fixed point of the spin-boson model. NRG 
indicates that the critical fixed point merges with the 
delocalized one as s — > + , thus we expect that the ex- 
pansion will allow access to the quantum phase transi- 
tion for small s. As shown below, the expansion is done 
about the lower- critical "dimension" s = 0; it yields an 
interacting fixed point, and mean-field critical behavior 
for small s does not obtain. For convenience we assume 
equal couplings, Aj = A, then the energy dependence of 
J{uj) is contained in the density of states of the oscillator 
modes u>i, and we have a a A 2 . 

How to set up a suitable RG expansion? Power count- 
ing about the free-spin fixed point, A — A = 0, gives 
the scaling dimensions dim [A] = (1 — s)/2, dim [A] = 1. 
Thus, both parameters are strongly relevant for small s. 
A better starting point is the delocalized fixed point, cor- 
responding to finite A. Eigenstates of the impurity are 
— > x ) and I <— x ), with an energy splitting of A. The low- 
energy Hilbert space contains the state | —> x ) only, and 
interaction processes with the bath arise in second-order 
perturbation theory, proportional to k = A 2 /A. Power 
counting w.r.t. the A = limit now gives dim[Ko] = — s, 
i.e., kq is marginal at s = 0, indicating that an e-type 
expansion for small s is possible. 

To lowest order, the RG can be performed within the 
low-energy sector, i.e., for the kq vertex (Fig. lb) and 
the propagator for the | — > x ) state. Consequently, our 
approach is valid as long as A <C A,ui c . We introduce a 
renormalized coupling k according to kq = A~ s k where 
A is the running cutoff, A = ui c initially. The one-loop 
beta function can be derived using the familiar momen- 
tum shell method, i.e., by successively eliminating high- 
energy bath bosons. To one-loop order only the coupling- 
constant renormalization in Fig. lc enters, and we obtain 

(3{k) = -SK + k 2 . (9) 

Besides the stable delocalized fixed point k = this flow 
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FIG. 1: Feynman diagrams occurring in the perturbation the- 
ory for the spin-boson model. Full/dashed lines denote the 
propagators of the | — * x ) and | <— x ) impurity states - the two 
states are separated by a gap of size A. The wiggly line is the 
local bath boson G\ oc - a) Interaction vertex A. b) Interaction 
vertex kq in the low-energy sector, c) One-loop renormaliza- 
tion of k. d) Diagrams for the local susceptibility xioc- 
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FIG. 2: Left: NRG data for the order parameter exponent 
1//3. Right: Magnetization exponent 1/8 from NRG, together 
with the analytical RG result (12) [solid line]. The dashed 
lines are the mean-field results f3 = 1/2, 5 = 3. 



equation displays an infrared unstable fixed point at 

k* = s + 0(s 2 ) (10) 

which controls the transition between the delocalized and 
localized phases. No (dangerously) irrelevant variables 
are present in this theory, so we conclude that the critical 
fixed point (10) is interacting [17]. 

We proceed with the calculation of critical exponents. 
Expanding the RG beta function around the fixed point 
(10) gives the correlation length exponent 

l/V = S + 0(s 2 ), (11) 

i.e., v diverges as s — ► + , as characteristic for a lower- 
critical dimension. Parenthetically, we note that the RG 
structure for s — > l~ is also similar to that near a lower- 
critical dimension: The line of second-order transitions 
for < s < 1 terminates in a Kosterlitz-Thouless transi- 
tion at s = 1 and is thus bounded by two lower-critical 
"dimensions" - a similar situation was recently found in 
the pseudogap Kondo problem, which, however, is in a 
different universality class [18]. As usual for RG expan- 
sions around a lower-critical dimension the present RG 
can only capture one of the two phases (the delocalized 
one), whereas run-away flow occurs on the localized side. 

The exponents associated with the local field e can 
be obtained in straightforward renormalized perturba- 
tion theory. To calculate observables, the diagrams arc 
written down using the original model with couplings A 
and A and both | — » x ), | <— x ) states. The perturbation 
theory turns out to be organized in powers of A 2 / A, as 
expected. Some of the relevant diagrams are displayed 
in Fig. Id, details will appear elsewhere. Restricting our- 
selves to the lowest-order results for the disordered and 
quantum critical regimes we find 

7 = 1 + C(s), x = y = s + 0{s 2 ) 1 
1/5 = l-2s + 0(s 2 ). (12) 

Interestingly, we are able to derive an exact result for the 
exponents x, y, employing an argument along the lines 



of Refs. 19, 20, based on the diagrammatic structure of 
Xioc- We obtain 

x = y = s . (13) 

(Notably, y = s was found to be the exact decay exponent 
of the critical spin correlations in the long-range Ising 
model for all s [9, 11].) Hyperscaling yields 5 = (1 + 
s)/(l — s), consistent with the lowest-order result (12). 

As an aside, we note that at s — the bath coupling 
is marginally relevant. Therefore the impurity is always 
localized as T — > 0, with a localization temperature given 
by T* = w c exp(-Aw c /A 2 ). 

Spin-boson model: Numerical results. We have per- 
formed extensive NRG calculations [8] to evaluate the 
critical exponents of the spin-boson transition, see 
Refs. 7, 8 for numerical details [21]. Results are shown in 
Figs. 2 and 3: these exponents obey hyperscaling includ- 
ing x = y (co/T scaling). They are in excellent agreement 
with the small-s RG expansion, but at variance with the 
exponents of the long-range Ising model: the mean-field 
predictions are /3 = 1/2, 5 = 3 which are clearly violated 
by our results in Fig. 2. 

Within error bars the mean-field exponents are real- 
ized at s = 1/2. Further, the one-loop results for v (11) 
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FIG. 3: Left: Numerical data for the correlation length ex- 
ponent \/v obtained from NRG, together with the RG result 
(11). Right: Susceptibility exponent x from NRG, together 
with the RG result (13). Here, the mean-field results coincide 
with the lowest-order perturbative ones. 
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and 7 (12) appear to be exact for all < s < 1/2, and 
logarithmic corrections to the power laws are observed 
at s = 1/2. This suggests that the spin-boson transi- 
tion does change its character at s = 1/2, but the critical 
fixed points for both < s < 1/2 and 1/2 < s < 1 are 
interacting and obey hyperscaling (the latter one being 
equivalent to that of the classical Ising model). 

Discussion. We have proven that the naive quantum- 
classical mapping fails for the sub-ohmic spin-boson 
model: Using a novel RG expansion around the delo- 
calized fixed point, we have shown that the quantum 
transition at0<s<l/2is controlled by an interact- 
ing fixed point, whereas the corresponding classical long- 
range Ising model shows mean-field behavior. Thus, the 
spin-boson problem for s < 1/2 is equivalent neither to 
the classical Ising model nor to the corresponding (quan- 
tum or classical) 0(1) </> 4 theory [10]. In physical terms 
the inequivalence can be traced back to the different dis- 
ordered (delocalized) fixed points in the two situations 
(expansions around these fixed points are suitable to ac- 
cess the critical behavior for small s): In the quantum 
model the transverse field fully polarizes the spin in x 
direction (which can be viewed as a "condensate" of spin 
flips), whereas the high-temperature limit of the classical 
Ising model is simply incoherently disordered. 

The inequivalence may come unexpected - so where 
does the quantum-classical mapping fail? Formal proofs 
of the mapping using transfer matrices rely on the short- 
ranged character of the interaction [1] . For general inter- 
actions a Trotter decomposition of the quantum partition 
function is employed where the imaginary axis of length 
/? = 1/T is divided into N slices of size At = (i/N, lead- 
ing to an Ising chain (8) with N sites. This procedure 
is exact when the limits At — > and /3 — ► oo are taken 
in this order. However, the limit At — > leads to a di- 
verging near-neighbor coupling in the term Hsr of the 
classical Ising model (8) [6]. This may in fact change the 
critical behavior of H c i (8), as a classical model with finite 
couplings arises upon taking /3 — > oo first. In other words, 
the quantum and classical problems are only equivalent 
if and only if the low-energy limit of H c \ is independent 
of the order of the two limits At — > and /3 — * oo [6] . As 
we have proven the inequivalence of TCsb and Ti. c i (with 
finite couplings) we conclude that the two limits cannot 
be interchanged for s < 1/2. 

A recent paper [22] investigated a SU(A^)-symmetric 
Bose-Fcrmi Kondo model in a certain large- N limit and 
found a critical fixed point with ui/T scaling for all s. The 
authors argued that the apparent failure of the quantum 
classical mapping is due to the quantum nature of the 
impurity spin. As discussed above, this alone is not suffi- 
cient: for short-range interactions in imaginary time the 
mapping can be proven to be asymptotically exact [1]; 
long-range interactions are essential. 

Our results suggest that some conclusions drawn in 
the past for effectively long-range Ising systems on the 



basis of the quantum-classical mapping have to be re- 
examined. Further, we envision that our novel e expan- 
sion will have applications for various quantum impu- 
rity problems, e.g., two- level systems coupled to multiple 
baths, Bose- Fermi Kondo models [20, 22], and also for 
quantum dissipative lattice models [23]. 

We thank D. Griineberg, H. Rieger, A. Rosch, S. 
Sachdev, Q. Si, M. Troyer, T. Vojta, and W. Zwerger for 
discussions. This research was supported by the DFG 
through SFB 484 (RB) and the CFN Karlsruhe (MV), 
and by the Alexander von Humboldt foundation (NT). 
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Erratum: Quantum phase transitions in the sub-ohmic spin-boson model: Failure of 
the quantum— classical mapping, Phys. Rev. Lett. 94, 070604 (2005) 

Matthias Vojta, Ning-Hua Tong, and Ralf Bulla 

We point out that the interpretation of data obtained by the Numerical Renormalization Group 
(NRG) method in Phys. Rev. Lett. 94, 070604 (2005) was incorrect: Two different effects precluded 
the reliable determination of critical exponents. An analytical argument was as well incomplete. 
We discuss consequences for the quantum phase transition scenario. 



Ref. 1 used results from NRG to argue that the critical 
properties of the spin-boson model with bath exponent 
s < 1/2 are not those of mean-field theory. The difference 
is in the order-parameter exponents (3, S and x. 

By now, we have convinced ourselves that the expo- 
nents in question could not be reliably extracted from 
the NRG calculations, if the ordered phase is controlled 
by an operator which is dangerously irrelevant at criti- 
cally. The two relevant issues will be discussed in turn. 

(A) Hilbert-space truncation: The bosonic NRG is 
presently unable to describe the physics of the localized 
fixed point for s < 1, due to the truncation of the bosonic 
Hilbert space [2] . From the absence of convergence prob- 
lems at the critical fixed point we nevertheless concluded 
[1] that all critical exponents can be extracted from NRG. 
While correct for an interacting fixed point, this assertion 
is erroneous for exponents [3 and 5 of a Gaussian fixed 
point, because, in the presence of a dangerously irrelevant 
variable, those exponents are not properties of the criti- 
cal fixed point and its vicinity, but instead properties of 
the flow towards the fixed point controlling the ordered 
phase. (This is what happens above the upper-critical 
dimension of a standard </> 4 theory.) 

(B) Mass flow: The NRG algorithm integrates out the 
bath iteratively. In the order-parameter language, this 
leads to a mass renormalization along the flow, arising 
from the bath (not from non-linear interactions), i.e., at 
any NRG step only a partial mass renormalization has 
been taken into account. As a result, the critical flow 
trajectory in NRG is that of a system with an additional 
mass set by the current NRG scale (i.e. by the tempera- 
ture T), and, effectively, the system is not located at the 
critical coupling for any finite T. For an interacting fixed 
point, the missing mass correction from the bath and the 
interaction-generated mass have the same scaling, hence 
no qualitative error is introduced. In contrast, above the 
upper-critical dimension, the interaction-generated mass 
vanishes faster as T — ► 0, and the bath mass correction 
dominates all observables calculated along the critical 
flow trajectory, like the exponent x. Note that this con- 
ceptual "mass-flow" problem of NRG arises for all baths 
without low-energy particle-hole symmetry. 



Consequently, the NRG-extracted exponents 0, 8, and 
x quoted in Ref. 1 were unreliable. The same applies to 
the analytical argument given in Ref. 1: The RG equa- 
tions (9-11) were correct, but the following analysis relied 
on the absence of dangerously irrelevant variables as well. 

Analyzing (A) in detail, we have found that careful 
NRG calculations performed with a large number 
of bosonic states in principle allow to extract the cor- 
rect exponents (3 and 5. While the truncation always 
dominates asymptotically close to criticality, the physi- 
cal power laws can be seen on intermediate scales. As an 
example, we show in Fig. 1 the non-linear field response 
at criticality for s = 0.4. For Nj, = 50 a crossover from 
5 « 2.4 to the mean-field value 5 = 3 upon increasing the 
field is visible, before the quantum critical regime is left. 

Taking into account recent numerical work based on 
Quantum Monte Carlo and exact diagonalization tech- 
niques [3, 4], which confirmed mean-field behavior in the 
sub-ohmic spin-boson model for s < 1/2, we conclude 
that the quantum-to-classical mapping is valid here. 
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FIG. 4: Magnetization curve of the critical spin-boson model 
for s = 0.4. NRG parameters are A = 4, (N s ,N b ) = (40, 12) 
and (60,50) [2]. The lines are power laws with exponents 1 
(dotted), 1/2.4 (dashed), and 1/3 (dash-dot). 



